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In a heavy fermion system, there exists the anomalous 
HaU effect caused by locahzed /-orbital freedom, in addition 
to the normal Hall effect due to the Lorentz force. In 1994, we 
found that the Hall coefficient caused by the anomalous Hall 
effect (-Rh'^^) is predominant and the relation oc 
(p is the electrical resistivity) holds at low temperatures in 
many compounds. In this work, we study the system where 
the magnetic susceptibility is highly anisotropic due to the 
strong crystalline electric field on /-orbitals. Interestingly, we 
find that is nearly isotropic in general. This tendency 

is frequently observed experimentally, which has casted suspi- 
cion that the anomalous Hall effect may be irrelevant in real 
materials. Our theory corresponds to corrections and gener- 
alizations of the pioneering work on ferromagnetic metals by 
Karplus and Luttinger. 
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In many heavy fermion (HF) systems, the Hall co- 
efficient (i?H) in a paramagnetic state shows universal 
temperature dependence, reflecting the crossover of the 
electronic property: As temperature increases, i?H also 
increases rapidly, and decreases after showing a peak at 
Tq. (we call Tq the 'coherent temperature'). The 

peak value i?H(2o) is apparently enhanced: it reaches 
^ 10^ times that of a normal metal. The origin of this 
interesting phenomena should be confirmed. 

Today, some theoretical explanations have been pro- 
posed in terms of the anomalous Hall effect (AHE), which 
is caused by the /-orbital freedom. In 1994, we 

proved the existence of the AHE in the periodic An- 
derson model (PAM) with orbital degeneracy, using the 
Kubo formula and the Fermi liquid theory. The ob- 

served i?H in HF systems is given by i?H — 
where is the anomalous Hall coefficient we found 

and /iji^f™ is the normal Hall coefficient caused by the 
Lorentz force, which is ternperature independent below 
the coherent temperature. [^,0 The AHE naturally ex- 
plains the sign and the magnitude of i?H in HF sys- 
tems, and it shows a simple relation at low temperatures, 
~ a ■ p^, where p is the observed electrical resis- 
tivity, (see Fig. 11 of ref. Q.) The mechanism of our 
AHE is essentially similar to that found in ferromagnetic 
metal by Karplus and Luttinger. [^HTl|] 

In this work, we extend this theory to an anisotropic 



system with a crystalline electric field (CEF). Figure |l| 
shows anisotropy of Rh and x in CeRuaSia |^], in which 
the relation i?H oc holds at low temperatures, as shown 
by Fig. ||. While x is extremely anisotropic, Rh is nearly 
isotropic. This appears quite paradoxical because AHE 
is the magnetic contribution and i?^^^ is enhanced by x 
in the spherical system M. We calculate the direction- 
dependence of both i?y and x under the tetragonal 
CEF, and find that i?^^^ is nearly isotropic even if x is 
highly anisotropic, in general. This study confirms that 
AHE is predominant at low temperatures in many HF 
compounds. 

The history of the study of AHE is very long, which 
started on the ferromagnetic metals. |jll|-jl^. Here, we 
comment on the other theory on AHE in paramagnetic 
HF materials given by Fert and Levy, Q which is the im- 
provement of ref. . They studied the impurity Ander- 
son model (J = 5/2) with an additional scattering poten- 
tial (1 = 2). They calculated i?^"^ using the Boltzmann 
equation for higher temperatures (T > Tq), where each 
/-electron site is almost incoherent. This mechanism is 
essentially the same as the theory of AHE in dirty fer- 
romagnetic metal, given by Smit. [ p^ They concluded 
that cx X ■ P for T > Tq. Here we comment that 

their theory on AHE are essentially different from ours, 
which cannot be derived by the Boltzmann equation. For 
example, the simple extrapolation of our theory predicts 
that (X X in the incoherent region (T > Tq), which 

is independent of p. Both of them may co-exist. But, it is 
definite that our AHE predominates at least for T < Tq. 

Here, we study Ce-compounds, where /^ configuration 
is realized. Due to strong l-s coupling, we have only to 
treat J = 5/2 sector, which has sixfold degeneracy. In a 
tetragonal crystal, this degeneracy splits into three levels 
by CEF, |r) = I ± 1), I ± 2), I ± 3), where | + F) and 
I — r) form a Kramers doublet. ( Now, is not a good 
quantum number.) In this work, we discuss the orbitally 
degenerate (J — 5/2) PAM under a tetragonal CEF as 
follows: 
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where cj^.^ iflr) is the creation operator of the conduc- 
tion electron (/-electron), and ek (-Epp/) is the spec- 
trum for conduction electrons (/-electrons). In a tetrag- 
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onal crystal, the /-orbital eigenstate |r) is given by 
|r) = Y.mOtm\J = ^I'^,Jz = M), where Otm is an 6 x 6 
orthogonal matrix and is given by 



0±l,±5/2 — 0±3^±3/2 — a, 



(2) 
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and other elements are zero. In a cubic crystal, a = 
•\/5/6 = 0.913. (In a hexagonal crystal, a — 1.) Then, 
we can easily verify that Vkro- — X^a/ '^rM VkMo- and 
-^rr' = Y.mm'OtmOt'M'Jmm' {l'^x,y,z), J'^ being 
the angular momentum along 7. (Note that VkMo- = Vq • 

\/4V3Em-'^\/(7/2-Ma)/75™,M-./2-i^i^3(^^fc,<Pfc) is 
derived for a spherical case, where Yi^^{9k,<fk) is the 

spherical harmonic function, jl^) Finally, ii'/p/ under 



the magnetic field H is given by Efj,, 
g^B '^-f J-pT' ■ -ff ^-order, where g is the Lande's 

(/-factor ( (7 = 6/7 for J — 5/2 ), and we put fiB — ^ here- 
after. For the time being, we limit our study to the case 
of /7 = 0, or to the mean-field (Gutzwillcr) treatment: 

Vkro- ^/z^ ■ VkFo- and sf ^ = zp ■ {eI - fi) + fi, 
where is the renormalization constant. [ p5[ The elec- 
tronic structure of orbitally degenerate J — 5/2 PAM is 
shown in Fig. |^, where i?k is the quasiparticle spectrum 
for k, given by the pole of Green functions. 
In the mean-field approximation, x is given by 
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where xp and xv are called Pauli susceptibility and Van- 
Vleck susceptibility, respectively. Both xp and xv take 
the positive value After a long calculation, we get 

the expressions in the case of H || 7 (7 = a;, z): 
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where p5^(e) represents the DOS for conduction electrons, 
—D is the bottom of E^, and 
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Q!kr(T(e) = Vl<;r(T/(e - -E^)- 
ak(e) = l + El"kr<T(e)P 
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Note that dE\s_/dek — ak(-Bk). The i?k-integration in (||) 
can be done easily when p^(e) is constant. In the case of 
no CEF (i?p = E^ for any F), the relation xp ~ Xv is 
derived by (^ and (||). In contrast, in the strong CEF 



case (£"/ ^ Ej,El), xp shows strong anisotropy and 
Xp 3> Xv is realized within the mean-field approxima- 
tion. In reality, it is not easy to investigate the accuracy 
of the above mean-field results. In the case of no CEF, 
this problem is studied by the Fermi liquid theory, \lt] 
or by the second-order perturbation theory w.r.t. U, |17] 
and the relation xp ~ Xv has been confirmed now. 

On the other hand, the Hall conductivity due to AHE 
is given by eq. (2.16) in Ref. Q, using the Kubo formula. 
While we take the effect of CEF into account, we neglect 
the vertex corrections for currents. Under the magnetic 
field H II 7 (■j = x,y, z), the anomalous Hall conductivity 
t^a™ (ea/37 = l) i^ given by 
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The derivative of the numerator in (^^ gives the anoma- 
lous transverse current. In (||), Ap (A^) is the damping 
rate of the /-electrons (conduction electrons) , given by 



Ap 
At 



-ImI]p(/^ + iO) 
El"kPo(M)r Ar, 
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where Ep(w) is the self-energy of /-electrons with F, 
which is k-independent within the local approximation. 
Here, Sr('^) is brought by the Coulomb interaction U. 
Within the second order perturbation and local approx- 
imation, taking the Pauli principle for F into account, it 
is easy to show that 
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for any T at low temperatures (T <C Tq), where p^ip) = 
Skcr l'^kr(T(M)P ■ PklA*) is the density of states for /- 
electrons with F at the Fermi energy. The relation (|l4| ) 
tells that Ar / Ap' is temperature- independent for any F 
and F' at low temperatures, which will also be expected 
in the coherent region [T <Tq). As a result, cr^|^^ given 
by is temperature-independent in the coherent region, 
where x remains nearly constant and Ar <IC {E^ — ii) is 
realized. Q The definition of the Hall coefficient for H || 7 
is i?H = {<yai3/H)/{aaaCFpi3) (whcrc ea/37 = l). Thus, the 
relation oc is confirmed, as shown in Fig. |^. 

The theory of Karplus and Luttinger contains a serious 
inconsistency, pointed out by ref. In ref. we 

solved it by taking the dissipation of quasiparticles into 
account. The new aspects of the AHE we find are as 
follows: (i) We derived the exact form of cr^p^^ ^^'^ 
found that it depends on the ratio Ar / Ap' , which differs 
greatly from 1 in the anisotropic system. This means 
that cr^p^ depends on the origin of the quasi-particle 

damping, (ii) In our theory, the orbital moment {gJ^) 
plays an important role, and the l-s coupling term is not 
essential. 

Following are the results of our numerical calculations 
on X and c^^^/ff, given respectively by and (||). 
Now we put zy — \. Here, we study the case A1-2, A2-3 > 
0, where A1-2 = - e( and A2-3 = E^ - e(. Taking 
the experimental data on CeRu2Si2 into account, we put 
A2-3 = 3- Ai-2 ||2^. For simplicity, we treat A1-2 and the 
parameter a independently. Other parameters we use are 



E 



f 



-0.2, 



0.23, and Vq = 0.2. We also put Ck = 
— H- 2 • (k/Tr)"^, where p^. — ■ dek/dk is k-independent. 
These selections (and the shape of the Fermi surface) is 
to be irrelevant to the conclusion qualitatively. Figure ^ 
shows the a-dependence of x^/x^ and o-^^^/c^y^^ in the 
case of A 1-2 = 0.3. We can recognize the isotropic feature 
of the anomalous Hall coefficient, which is quantitatively 
affected by the form of the damping rate. (When A1-2 = 
A2-3 = 0, both X and (Tab are isotropic and independent 



of a.) Figure || shows the dependence of 



X 



^AHE 



X 5 ^ xy 

and cTy^^ on A1-2, in the case of a = 0.96. We can see 
that -Rh is less affected by A1-2 than x- As a result, 
-Rg^^ can take a large value even in strongly anisotropic 
systems. Especially, tends to be isotropic for larger 

Ai-2 (Ai-2» \E(~fi\). 

Such a difference of the anisotropy between on x and on 
cr^p^ can be understood naturally as follows: xp: given 
by (^, comes mainly from the lowest Kramers doublet 
(say, F — ±1). On the contrary, higher Kramers doublets 
are indispensable for a^p^ , given by (^) , because AHE is 
the multi-band contribution. Thus, cr^^^ is averaged 
and is nearly independent of directions. In conclusion, 
we find that the anisotropics of x and R^^^ can be quite 
different. In the strongly anisotropic system due to the 



(tetragonal) CEF, in general, x is strongly anisotropic 
and R^^ tends to be isotropic. 

Finally, we point out some remaining problems. The 
present model may be too simplified for quantitative 
study in that it contains only one spherical conduction 
band ek. Besides, some kinds of vertex corrections (or 
1 /zr) omitted in this work should be studied for complete 
analysis. Moreover, the consideration on the Umklapp 
processes will be necessary for some compounds: Refer- 
ences show that axx oc ■ Cum and a^™/H oc 
{Cvm}^, thus i?H (X Al ■ {Cum}°, where the coefficient 
Cum reflects the Umklapp processes for electron-electron 
scattering, which dominate the various transport phe- 
nomena at low temperatures. |0 (Note that the Umk- 
lapp processes are irrelevant to i?H-) Cum ~ 0(1) is 
satisfied in the usual case. But Cum ~ 10 will be re- 
alized in CeRu2Si2, which shows large a^^^/H at low 
temperatures, (see Fig. || or ref. [Q.) At last, we com- 
ment on the existence of the AHE in other kinds of met- 
als. For example, at low temperatures, i?H in Sr2Ru04 
shows the temperature dependence similar to that in HF 
compounds, jl^ We find that the relation i?H oc holds 
for T < 6K. pO| We believe that the AHE discussed in 
this work is general to many non-magnetic metals with 
orbital- freedom. 
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FIG. 3. EflFcctive band structure for J = 5/2 PAM under 
the tetragonal CEF. n is the Fermi energy. 



20.0 



10.0 



0.0^ 



0.92 



0.94 



0.96 



0.98 



FIG. 4. a-dependence of the anisotropy of x and a^^^ . {a 
characterizes the CEF introduced by eq. (2).) Here, we set 
Ai-2 = 0.3 and \E( — = 0.03. In comparison, we also 
^ /(^x™ for the case of Ar = constant, instead of 



show ( 

eq. (14), by staxs (dash-dotted line). 




FIG. 5. AcEF-dependence of x^, x'^ , o-^™ (H || z) and 
™ (H II x). Here, we set a = 0.96 and \E( - m| = 0.03. 

FIG. 2. Rh vs for CeRu2Si2. This fi gure is borrowed 
from ref. [3]. 
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